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Abstract 

Ground state properties of trapped Bose condensate with repulsive interaction in non-homoge- 
neous gravitational field are studied. Spatial structure of Bose condensate and its momentum 
distributions in 3-D anisotropic trap are considered by the solution of the modified non-linear 
Schrddinger equation. The results are compared with the corresponding properties of condensate 
in a harmonic trap without gravitational field. 

PACS numbers: 03.75.Fi, 05.30. Jp, 32.80.Pj 

Bose-Einstein condensation (BEC) is an outstanding example of a case where the properties of 
quantum matter can be studied in laboratory experiments. Since the theoretical discovery of BEC 
in the mid-1920s, most studies in this area were concentrated on the properties of homogeneous 
condensates with effectively infinite spatial extents. However, experimental success in achieving 
and observation of Bose condensation in trapped dilute atomic gases, which was done during the 
last five years, provided the possibility for study of the properties of an inhomogeneous confined 
condensate A trapped Bose gas exhibits features that are impossible to observe in a 

homogeneous medium. One of the remarkable features of these experiments is the observation of 
Bose condensation in coordinate as well as in momentum space. The localization in space of the 
condensed fraction of a low temperature Bose gas allows the detection of the condensate and to 
estimate the size of the BEC sample, to find its shape and momentum distribution of its atoms 
Moreover, the ability to manipulate such coherent matter provided additional stimuli for 
experimental and theoretical study of this phenomenon. The properties and behavior of a BEC 
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in a trap depend on the presence of various factors. One of these factors is the interaction of 
the trapped atoms with a gravitational field. Gravity has a dominant influence on the evolution 
of atoms released from the trap. The dominance of gravity over the thermal effects has been 
successfully used in atom cloud gravimeters 0-0- The use of BEG as a gravity and gravity 
gradient sensor opens wide perspectives for the creation of highly sensitive gravimeters and 
gravity gradiometers [^]. In this paper we present quantitative results on the influence of a non- 
homogeneous gravitational field on the behavior of a Bose condensate with repulsive interactions 
in a 3-D anisotropic trap. We determine the properties of a BEG in a gravitational field by 
studying its geometry, momentum distribution, and dynamical properties. We concentrate on 
the results obtained by the analysis of the solutions to the non-linear Schrodinger equation for 
trapped atoms in external gravitational field. Zero temperature solutions for the condensate are 
found using the Gross-Pitaevski (GP) equation Out objective is to provide a general 



framework for the analysis of the behavior of a confined Bose condensate which interacts with a 
gravitational field through the modified GP-equation and to investigate its properties. 

A BEG created in a trap is under the influence of different forces, and the properties of 
the condensate depend on the combination of three potentials. A magnetic trap generates a 
time-averaged harmonic potential, while the second effective potential is formed from the binary 
collisions of the atoms. Finally, the condensate is under the influence of an external gravitational 
field which gives rise to a gravitational potential. The properties of the condensate at zero tem- 
perature are described by GP theory [|13|-|15[ which in its standard form includes two potentials. 



Our aim is to modify the GP equation in order to take into account gravitational interactions. 
The many body Hamiltonian describing N interacting bosons with atomic mass m is given by 
the equation 



H = j rfx^+(x) 



-A + v;^t(x) 



^(x) 



+^ j dx'dxV^+(x')^/'+(x)y(x' - x)^/'(x')V'(x), (1) 
where '?/'(x) and are the boson field operators which annihilate and create a particle at 



the position x, respectively, and l^(x' — x) is two body interatomic potential. Time evolution 
for the field operator ip{'x,t) is found from the Heisenberg equation 



-—A + Kxt(x) 
2m 



dx.'ip^{x.' ,t)V{x.' - x)^/'(x', t) 



(2) 



Only binary collisions at low energy are relevant for the dilute gas at low temperatures. These 
collisions are well described by s-wave scattering processes and allow us to replace V{'k — x) 
with an effective interaction 

\/(x'-x) = K5(x'-x), (3) 
where the coupling constant Vo is related to s-wave scattering length a through 



Vo 



m 



(4) 



To overcome the problem of solving the full many-body Schrodinger equation for interacting 
systems we can use the mean field approach [0,[|T^. The mean field approach is based on 
separation of the condensate fraction of the bosonic field operator. For a non-uniform and time- 
dependent configuration we have 



7/'(x,t) = 0(x,t)+f (x,t), 



(5) 



where 0(x, t) is a complex function defined as the expectation value for the field operator 
(f){x,t)=< V^(x, t) >. As a result, the equation for the condensate wave function 0(x, t) is 
found using (0) and (0), and is written in the form 



_iLA + K.t(x) + K|0(x,t)P 
2m 



(x,t). 



(6) 



This non-linear Schrodinger type equation is called the Gross-Pitaevski (GP) equation. The 
validity of this equation is based on the conditions that the s-wave scattering length is smaller 
than the average distance between atoms, and that the number of atoms in the condensate 



is large. To obtain the GP-equation for the ground state, one can write the condensate wave 
function as 

0(x,t)=0(x)e-^'^*/\ (7) 

where the wave function 0(x) is normahzed to the total number of particles / (ix|0(x)p = A^, 
and /i is a chemical potential. Inserting into (|^), we find the equation for the wave function 

-^A + VU^) + 1!^|0(x)p1 0(x) = /i0(x). (8) 

The external potential for the confined atoms is the sum of the trap potential and a contribution 
due to the interaction with gravitational field: V^xt = Vtrap + Vg- The magnetic trap generates a 
harmonic confining potential of the form 

Vtrap{^) = y (c^^X^ + + ulz") , (9) 

where uji are the projections of the angular trap frequencies on the coordinate axis, and vector 
X = {x,y,z) defines the displacement from the trap center Xq = {xo,yo,Zo)- The gravitational 
interaction is described by ^^(X) = m$(X), where the potential for the given mass distribution 
is found as $(X) = —G J2i ^i/|Xj — X|. The potential of a gravitational field is an unobservable 
value. The intrinsic characteristics that uniquely characterize the gravitational field in Newtonian 
gravity are the first and the second spatial derivatives of the potential. They form the vector of 
gravitational acceleration Qi = —d^/dxi, and the gravity gradient (or Eotvos) tensor 

Each component of this tensor defines the acceleration difference along the direction j per unit 
separation along the direction k. The gravity gradient tensor is symmetric, and its trace is 
related to the local mass density by Poisson's equation. These conditions define the number of 
its independent components. In the local coordinate basis of the trap the gravitational potential 
can be written with the help of the vector gravitational acceleration and the gravity gradient 
tensor as 

<l>(xo, x) = $(xo) - gj{^o)x^ + ^TjkMx^x''. (11) 



We will locate the coordinate system of the trap in such a way as to obtain the components of 
the vector of gravitational acceleration as 



9x = 0, gy = 0, gz = -\g\, 



(12) 



In this approximation, the leading components of the gravity gradient tensor are {T^x, ^yy, Tzz}- 
The components of the gravity gradient are not independent. They are connected with each 
other by the Laplace equation which leads to T^^ + Tyy = — F^^. In the given coordinate system, 
the equation for the external potential is 



(ul + Txx) + {ul + Tyy^ + (uj + r^^) z'^ + mgz + m$o 



(13) 



Consider an axially symmetric confined potential with angular frequencies = uoy = uj±_ and 
ujz = Xuj±, where A is an asymmetry parameter. Then, the symmetry of the trap allows us to 
write down GP equation as 



where 



A + ^ 

2m 2 



+mgz + 



4:7ih a 



m 



(/ig -m<l>o)</)(x). 



(14) 



7^(r)= 1 + ^ 



2 I ' 
1 + ^\ 



7.^(r) = AMl + ^ . 



(15) 



In the equation ([TsD, the factors {•y^iT) , ■jyiT) , •jzi^)} , are functions of the diagonal components 



of gravity gradient tensor. In order to simplify the analysis of the equation ([l^ , we will scale 
the coordinates using a spatial scale parameter ( = {Ti/mui^y^'^ : 



xi = C"'x. 



(16) 
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We also introduce the scaled wave function 



(17) 



which satisfies the normalization condition 



rf'xi|0i(xi)P 



1. 



In the new variables, the dimensionless GP-equation (|T^) has the form 



-Ai + 7^(r)x? + ^liT)yf + ^liT)zf + 2Dig)z, 



+Mi|0i(xi,?/i,^i)p 0i(xi,?/i,2;i) = 2/ii0i(xi,yi,Zi) 



where 



Ui 



is the dimensionless interaction strength, 



-{fig - m$o) 



(19) 



(20) 



(21) 



m 



(22) 



is a dimensionless eigenvalue of the equation, and 

fiuj± ^ 

is a linear function of gravity. The dimensionless form of GP equation (^^ is suitable for 
integration by means of appropriate methods. The non-linear term in this equation does not 
allow us to obtain analytical solutions, and numerical computer simulation has to be used instead. 
However, it is possible to find analytical solutions in two limiting cases: non-interacting particles 
and strongly repulsive interactions. 
(a) Condensate without interaction 

If we neglect two-body collisions, solving the GP equation can be reduced to finding the solution 
of (|TD with ui=0 

Dig) ' 



(j)i{xi,yi,zi] 



2/ii0i(a;i, yi, zi) 



(23) 



where 



(24) 



The dimensionless ground state solution is 



zi) = Aexp 



exp 



2 

D{9)\ 



2 

7.(r) 



X 



Zl + 



2 \- ' 7Kr) 

where the normahzation constant is found from condition (|TBp as A = {'Jxlylz/'^^Y^^, and 

/ii(^7,r) = (7.(r) + 7,(r)+7.(r))/2. 

The ratio of the mean values of momenta can be found using the wave function 



(25) 



(26) 



MPx,Py,Pz) = exp 



2n' 



(27) 



From this it follows that the ratios of the mean momenta are 



<Pi> 



7.(r) 



\<PI> \lA^y \<Pl> \lyiX) 



<pi> 



7.(r) 



<pI> 



7x.(r) 



(2^ 



As consequence, the anisotropy of the velocity distribution depends on the gravity gradient. As 
the Bose condensate is characterized by the absence of thermal excitations, equations ( PBD reflect 
the result of atomic motion in the trapping potential only. The transverse and vertical widths 
of the Gaussian distribution ( [25| ) are 

1 o 1 o 1 



< xj >-- 



< vi >-- 



< zl > -- 



(29) 



27.(r)' 27,(r)' 27,(r)' 

and depend on the gravity gradient. The average < zi >= —D{g)/'~f'^{r) defines the shift of the 
condensate cloud along the 2;-axis. 

(b) Condensate in a repulsive limit 

The parameter ui of equation ([l9|) is proportional to the total number of atoms in the condensate. 



In the limit of large one can neglect the kinetic energy term (the so-called Thomas-Fermi 
approximation) and rewrite this equation as 



+ui\(j)i{xi,yzi)\^ (l)i{xi,yzi) = 2fii(f)i{xi,yi, zi). 



(30) 



The ground state solution of equation ( pOD is 



h[xi,yi,zi} 



1 

Ul 



2/2i - ^l{T)xl - ^'^{T)yl - 7^(1) (^z^ + 



Djgy 
7.^(r), 



(31) 



where Hi is found from the normalization condition (KW) as 



2/ii(^7,r) 



15 

15aN 



c 



(7.(r)7,(r)7.(r))ni 
(7.(r)7,(r)7.(r)) 



1 2/5 

2/5 



(32) 
the 



The chemical potential depends on the trapping frequencies entering the potential Vtrap, 
number of particles A^, and the components of the gravity gradient. The right hand side of 
equation (|3l|) defines the area of the condensate location. The boundary of the condensate cloud 
has the form of an ellipsoid with semiaxes 

R „ , R „ , R 



7x(r)' 



7.(r)' 



Rzig,T) 



7.(r)' 



(33) 



where R = y2jli{g, T). As follows from equation (pi]), the center of the ellipsoid is shifted along 
the z-axis by —D{g)/j'^(T). 

In the conclussion we will review the obtained results. The expression for wave function 
( p5D in the the approximation of weak interaction between atoms allows us to find statistical 
characteristics of the density distribution (such as first and second cumulants, moments, mean 



and variance along each axis). As follows from (^5]) there is a normal distribution along the 
z-axis as well as along the x and y-axis. The distribution along the z-axis is characterized by the 
moment mi = —D{g)/'y'^{r), and dispersion cr^ = (272(r))^^. Therefore, the dispersion o"^ varies 
with respect to the Tzz component of the gravity gradient tensor, and the density distribution 
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is broadened in the z direction due to tidal forces. The condensate is shifted in proportion to 
the gravitational acceleration by < >= nii = —D{g)/'-j1{T). Similarly, the dispersions of the 
distributions in x and y are given by cr^ = (27^.(r))~^ and cr^ = (27y(r))^^. This shows that the 
density distribution is squeezed in the transverse direction due to the gravitational interaction of 
the atoms with the gravitational field. The central density of this distribution is n(0) = |0(O)p = 
{N/C^"^) [{'^xij')ly{^)lzij'))/'^'^^^'^ and depends on the number of particles in the condensate and 
the components of the gravity gradient through the 7-factors. The anisotropy in the velocity 



distributions \j<vl> / <vl> = ^j,{T)/-f^{T) and ^<v^^> / <v^> = ^-f,(T)/-fy(T) for 
non- interacting condensates is given by the ratio of 7-factors (Eq. ^8]), and therefore is connected 
with the variation of the diagonal components of the gravity gradient. In spite of the axial 



symmetry of the trap potential, there is an asymmetry in the ratio J< Vy > / <v1>, which is 



given by the ratio of 7-factors y7y(r)/7a,(r). As follows from equation (|3TD , the shape of the 
condensate in the strongly repulsive limit also depends on the diagonal components of the gravity 
gradient tensor. In the strongly repulsive limit, the 3-D shape of the condensate is an ellipsoid 
with semiaxes given by equation (^5]). The ellipsoid is squeezed in the transverse (x — y) plane, 
and extended in the axial {z) direction due to the contributions of the gravity gradient. The 
ellipsoid is also shifted in the z direction due to gravity. The energy per particle can be found 
from the equation jl = dE/dN and is given by E/N = 2/7 jl, where the dimensional chemical 
potential is jl = [15aA^(7^(^)7y(^)7^(^))/C]^^^ The contributions to the energy E in 

this approximation are related to the interaction between particles, the gravitational interaction 
and oscillator energies. The density of particles n{x,y,z) = \(j){x,y, z)\'^ is found directly from 
equation (0). As follows from equation (|3T|) , the density distribution vanishes in classical turning 
points on the surface of ellipsoid. The central density is n{0) = 2fliN/ {(^ui), and depends on the 
number of particles and the gravity gradient. The size of the ellipsoid depends on the number 
of particles A^ (the semiaxis scale as N^^^) and the components of the gravity gradient. The 
components of the gravity vector and the gravity gradient tensor are defined by the source of 
the gravitational field. For example, for Earth gravity {g = 9.87 ms~^ and T^z = 3 x 10^^ -s"^). 
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the correction for 7z(r) depends on the frequency (Eq.|l5D and can be numerically estimated 
as T.Jul ~ 10"^ for uj, ~ 2QHz. 
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